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NOMENCLATURE 


A 

a COL 

Alat 

a long 

A n _i,a n ,A n+ i 

Ar 

Ai,A2,A% 

A \ , A 2 , A 3 

Ao,Ai,A2,A3,A4 

B 

B\,B2,Bz 

B\,B2,B3,B4 

B\,B2, B3, B4, Bq 

c 

Cxi Cy, C Z 

Co,C]_,C2 

Cl,C2,C3 

Ci,c 2 ,c s 

C\,C2,- " ,C\4 

V 

D 

D 

Bxi Dy, D z 


general coefficient for the sin(cjt) term 

least upper bound (l.u.b) for the Nonrotating System Collective Fourier 
coefficient magitude constraint function 
least upper bound (l.u.b.) for the Nonrotating System Lateral Cyclic Fourier 
coefficient magitude constraint function 
least upper bound (l.u.b.) for the Nonrotating System Longitudinal Cyclic 
Fourier coefficient magnitude constraint function 
least upper bounds (l.u.b.s) for the Rotating System JV-1, N, and N+ 1 
per rev Fourier coefficient magnitude constraint functions, respectively, 
defined in section 2.3 

general least upper bound (l.u.b) for a rate control magnitude constraint 
function 

least upper bounds (l.u.b.s) for the first, second, and third actuator vertical 
excursions, respectively, defined in section 2.3 
coefficients appearing in equation (14) as defined in equations (15) of 
section 2.4.2 

general coefficients of the first product term appearing in equation ( 53 ) of 
section 2 . 6.2 

general coefficient for the cos (ut) term 

coefficients appearing in equation (14) as defined in equations (15) of 
section 2.4.2 

general coefficients of the second product term appearing in equation ( 53 ) of 
section 2 . 6.2 

general coefficients appearing in equation (48) of section 2.6.1 
general coefficient for the sin(atf + 7 ) term 

cofactors of the {x — 37), (y — y\), and (z — z\) elements in the V matrix, 
respectively 

general coefficient appearing in equation ( 49 ) of section 2 . 6.1 
amplitudes of the first, second, and third actuator vertical excursions, 
respectively, defined in section 2.3 

coefficients appearing in equation (14) as defined in equations (15) of 
section 2.4.2 

general coefficients appearing in equation ( 55 ) of section 2 . 6.2 
matrix associated with the swashplate plane equation defined by 
equation (2) of section 2.4.1 

determinant of the % matrix defined by equation ( 6 ) in section 2 . 4.1 
the constant of proportionality for the kinetic energy T appearing in 
equation (43) in section 2.6 

constants appearing in the Canonical form of the swashplate plane equation 
[eq. (5) of sec. 2.4.1] which equal D/C X ,D/C y , and D/C z , 
respectively where D = DetiJV) 


IV 



Dq,D\,D 2 

d 0 ,d u d 2 ,'--,d 12 

E 

E\ ,E 2 ,E 2 

F 

F 

F s 


F\,F 2 ,- ■ ■ ,f l 


G 

g.l.b. 

G\, G 2 , • • • , G m 
H 


H c 

HHC 

H s 

i 

3 

K 


k 

Ke 


K p 

l 

l 

l.u.b. 

m 

m 

N 

O 


general coefficients appearing in equation (50) of section 2.6.1 
general coefficients appearing in equations (58) and (59) of section 2.6.2 
total mechanical energy defined by equation (38) in section 2.6 
amplitudes of the first, second, and third actuator vertical excursions, 
respectively, appearing in equation (17) of section 2.4.2.1 
general force unit 

constant appearing in equation (42) of section 2.6 which equals 
one-half of the effective spring constant of the mechanical system 
force acting on the mechanical system due to the effective spring constant 
appearing in equation (39) of section 2.6 
general coefficients of the sin(0 + (f>s k ) terms appearing in 
equations (A-l) and (A-10) of appendix A 
constant appearing in equation (42) of section 2.6 which is the reaction 
to the spring constant at zq 
greatest lower bound 

general coefficients of the cos(0 + <f) Ck ) terms appearing in 
equation (A-10) of appendix A 

constant appearing in equation (42) of section 2.6 which represents a 
reference energy equal to K{zqzqq — 2 Zqq) 
constant defined by equation (A-4) or (A-13) of appendix A which equals 
the sum of the phase angle cosine terms 
higher harmonic control 

constant defined by equation (A-4) or (A-13) of appendix A which equals 
the sum of the phase angle sine terms 
blade number; i = 1, 2, . . . , N 

index number for the Do, D\,D 2 ,-" , D\ 2 general amplitudes appearing 
in equations (58) and (59) of section 2.6.2 
effective spring constant of the mechanical system with units of 
F/l = m/t 2 appearing in equation (39) of section 2.6 
actuator number; k = 1,2,3 

constant of proportionality that relates the mechanical power (Pm) t0 the 
time rate-of-change of the total mechanical energy (E) appearing in 
equation (37) of section 2.6 

constant of proportionality that relates the hydraulic power (Pfj) to the 
mechanical power (Pm) appearing in equation (36) of section 2.6 
general length unit 

pitch link connecting arm length appearing in equation (29) of 
section 2.4.2.3 
least upper bound 

number of oscillations per revolution in section 2.2 
general mass unit 
number of blades 

origin of the X — Y — Z coordinate system defined in section 2.3 
(see fig. 1) 


v 



p 

Pl,P2,P3 

P H 

p m 

n 

R 


sc(-) 

T 

t 

t 

U 

V 

Vi,V 2 


w 
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X,X\,X2,X3 

Y 

y 

y,yi,y2,y3 

yo 

z 

2,21,22,23 

20 

20 

200 

200 

2Qi , 202 , 2O3 


general point on the swashplate plane (see fig. 1) 

points of intersection of the swashplate plane with the vertical axis of the first, 
second, or third actuator, respectively (see fig. 1) 
hydraulic power requirement defined by equation (36) of section 2.6 
generated mechanical power defined by equation (37) of section 2.6 
matrix associated with the canonical swashplate plane equation [eq. (5) 
of sec. 2.4.1] 

radius vector from the origin (O) to a general point (P) on the swashplate 
plane (see fig. 1) 

radius magnitude of the right circular cylinder described by the rotating 
pitch links (see fig. 2) 

composite trigonometric function which denotes either sin(-) or cos(-) 
kinetic energy of the mechnical system appearing in equation (38) of 
section 2.6 
general time unit 
time 

control 6-vector [E\ E 2 E% (f>\ 02 0 3]^ defined by equation (18) of 
section 2.4.2.1 

potential energy of the mechanical system appearing in equation (38) of 
section 2.6 

potential energy of the mechanical system at states 1 or 2, respectively 
appearing in equation (39) of section 2.6 
general sum of monofrequency sinusoidal terms defined by 
equation (A-l) or (A-ll) of appendix A 
first coordinate axis of the X — Y — Z coordinate system defined in 
section 2.3 (see fig. 1) 

component of point P, P\,P 2 , or P3 along the X-axis (see fig. 1) 
second coordinate axis of the X — Y — Z coordinate system defined in 
section 2.3 (see fig. 1) 

general sum of monofrequency sinusoidal terms of defined by equation (30) 
of section 2.5.1 and by the equation (33) of section 2.5.2 
component of point P, Pi,P 2 , or P3 along the Y-axis (see fig. 1) 
average value of the sinusoidal y defined in section 2.5 
third coordinate axis of the X — Y — Z coordinate system defined in 
section 2.3 (see fig. 1) 

component of point P,Pi,P 2 , or P3 along the Z-axis (see fig. 1) 
reference actuator length defined in sections 2.3 and 2.6.1 
vertical displacement (z) at which the effective spring force is zero in 
equation (4) of section 2.6 

reference value of z which corresponds to the reference value of 6 , 6q 
in equation (29) of section 2.4.2.3 

reference value of z at which V\ is defined in equation (39) of section 2.6 
average value of the sinusoidal z\, z 2 , or 23, respectively defined 
in sections 2.3 and 2.6.1 
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0 

P 
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7 

Cl, <2 
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e 

e 

6c 

Qcol 

°LAT 

Qlong 

On-I^N^N+1 

e s 
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c 
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0S 
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01 , 02,03 


0 

01,02 

U) 

n 

Subscripts 

c 

COL 

E 

H 

i 

3 
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general angular argument 
general phase angle 

phase angle equivalent of the H s and if c coefficients appearing in 
equation (A-5) or (A-12) of appendix A 
general phase angle 

angular argument equal to either ( N — l)f)£, NQt, or (N + l)f It in 
section 2.1 

phase angle equivalent of the A and B coefficients appearing in equation (30) 
of section 2.5.1 

phase angles appearing in equation (49) of section 2.6.1 
general Fourier coefficients 

blade pitch angle defined in section 2.4.2.3 (see fig. 3) 
general angular argument appearing in equations (A-l) and (A-10) of 
appendix A 

general Fourier coefficients for cosine terms in section 2.1 
general Fourier coefficients for the collective terms in section 2.2 
general Fourier coefficients for the lateral cyclic terms in section 2.2 
general Fourier coefficients for the longitudinal cyclic terms in section 2.2 
general Fourier coefficients for the N-l, N, N+l harmonic terms in 
section 2.1 

general Fourier coefficients for sine terms in section 2.1 
reference value of the blade pitch angle defined in section 2.4.2.3 
(see fig. 3) 

dummy integration variable corresponding to z in equation (40) of section 2.6 
phase angles of the cosine terms of in equation (A-10) of appendix A 
phase angle used to define the relation between the rotating system and the 
nonrotating system in equation (19) of section 2.4.2.2 
phase angles of the sine terms in equation (A-10) of appendix A 
phase angle of the sinusoidal motion of the first, second, and third actuators, 
respectively in equations (16) and (20) of sections 2.4.2.1 and 2.4.2.2, 
respectively 

phase angle equivalent to the 6 S and 6 C coefficients appearing in the equation 
defining y in section 2.1 

phase angles appearing in equation (50) of section 2.6.1 
general frequency coefficient of t in sinusoidal terms 
rotor rotation rate 

cosine terms 

collective terms in section 2.2 
energy in section 2.6 
hydraulic processes in section 2.6 
blade number; i = 1, 2, . . . , N 

index for the Dj coefficients appearing in equations (58) and (59) of 
section 2.6.2 

actuator number; k = 1, 2, 3 
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LAT 

LONG 

M 

M 
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N 

P 

s 

s 

X 

y 

z 
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0 

00 

0,1,2 


0 , 1 , 2, -..,12 

1,2 

1,2 


1,2,3 
1 , 2 , 3 , 4, 5 

1 , 2 , 3, 14 


number of monofrequency sine terms in equation (A- 10) of appendix A 

lateral cyclic terms in section 2.2 

longitudinal cyclic terms in section 2.2 

mechanical power and mechanical processes in section 2.6 

number of monofrequency cosine terms in equation (A- 10) of appendix A 

least upper bound (l.u.b.) for a constraint function 
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power in section 2.6 

sine terms 

rotating system, pitch links, blades, etc., in section 2.4.2.2 
(x-xi) element in matrix V in section 2.4.1 
(y-yi) element in matrix V in section 2.4.1 
(z-zi) element in matrix V in section 2.4.1 
reference actuator length 

vertical displacement of which the effective spring force is zero in section 2.6 
reference value of 2 corresponding to the reference value of 6, 6q 
index values for the Co, C\, and C2 coefficients appearing in equation (49) of 
section 2 . 6 . 1 , and for the Dq, D\, and D2 coefficients appearing in 
equation (50) of section 2.6.1 

values for the j index of the Dj coefficients appearing in equations (58) and (59) 
of section 2 . 6.2 

states at which the potential energy of the mechanical system is defined 
index values for the £i and £2 phase angles appearing in equation (49) of 
section 2.6.1, and for the and \&2 phase angles appearing in equation (50) 
of section 2 . 6.1 
index values for the actuators 

index values for the B\, B2, £ 3 , B4, and B 5 coefficients appearing in 
equation (48) of section 2.6.1 

index values for the C\, C2, C 3 , • • • , C 14 coefficients appearing in 
equation (55) of section 2.6.2 


SUMMARY 


An important class of techniques to reduce helicopter vibration is based on using a Higher Harmonic 
controller to optimally define the Higher Harmonic blade pitch. These techniques typically require solu- 
tion of a general optimization problem requiring the determination of a control vector which minimizes 
a performance index where functions of the control vector are subject to inequality constraints. Six 
possible constraint functions associated with swashplate blade pitch control were identified and defined. 
These functions constrain: 1) blade pitch Fourier Coefficients expressed in the Rotating System, 2) blade 
pitch Fourier Coefficients expressed in the Nonrotating System, 3) stroke of the individual actuators 
expresed in the Nonrotating System, 4) blade pitch expressed as a function of blade azimuth and actuator 
stroke, 5) time rate-of-change of the aforementioned parameters, and 6) required actuator power. The 
aforementioned constraints and the associated kinematics of swashplate blade pitch control by means of 
the strokes of the individual actuators are documented herein. 



1.0 INTRODUCTION 


An important class of techniques to reduce helicopter vibration is based on using a Higher Har- 
monic controller to optimally define the Higher Harmonic blade pitch. These techniques typically require 
solution of a general optimization problem requiring the determination of a control vector which min- 
imizes a performance index where functions of the control vector are subject to inequality constraints 
(see sec. 2.3.1 of ref. 1). When solving these problems, it is extremely important to remember that 
the “truthfulness” or “correctness” of the solution depends not only on the validity of the mathematical 
methods employed to obtain the solution, but also on the “truthfulness” or “correctness” of the problem 
definition. In deductive mathematics, as in the case classical Aristotelian deductive logic, it is necessary 
that BOTH the syllogistic deductive reasoning be VALID and the assumed propositions be TRUE in 
order to assert that the conclusion is TRUE. 

The “truthfulness of the propositions” is often overlooked or ignored when attempting to solve 
difficult mathematical optimization problems with sometimes unfortunate results. Obviously, approxi- 
mations to “truthful” propositions can be required in order to reduce the problem to a tractable form, 
but oftentimes these approximations are made because of misunderstandings of what the problem really 
is (i.e., the real definition of the performance index and constraints). Concern about the correctness 
of the constraint formulation for the helicopter vibration minimization problem prompted the analysis 
documented herein. 


2.0 TECHNICAL 

Six possible constraint functions associated with swashplate blade pitch control were identified and 
defined. These functions constrain: 1) blade pitch Fourier Coefficients expressed in the Rotating System, 
2) blade pitch Fourier Coefficients expressed in the Nonrotating System, 3) stroke of the individual 
actuators expressed in the Nonrotating System, 4) blade pitch expressed as a function of blade azimuth 
and actuator stroke, 5) time rate-of-change of the aforementioned parameters, and 6) required actuator 
power. The aforementioned constraints and the kinematics of swashplate blade pitch control referenced 
to the individual strokes of the actuators in the nonrotating system are described in the following 
subsections. 


2.1 Rotating System Fourier Coefficients 

If given the Fourier Coefficients for an entity expressed in the Rotating System, say fyv_i) c , 
0(N- l) s > 6n c i 0n s > 0(N+l ) c » °{N+i) s 

where 

6 is a Fourier Coefficient 
N C is the number of blades 
c refers to a cosine term 
s refers to a sine term 

consider a pair of monofrequency coefficients 8 C , 6 S 


2 


Let 


y = 9 S sin 7 + 0 C cos 7 


where 

7 is either (iV-l)Qi, or iVf2L or (N+\)Clt 
Q is the rotor rotation rate 

t is the time 


Noting that 


sin(o! ± 0) = sin a cos 0 ± cos a sin 0 
cos(o: ± 0) = cos a cos 0 ^ sin a sin 0 


Then 


or 


y = 


9 S \ . , ( 9 C 

sin 7 + 


where ip is phase angle 
and 

Then 

So 


\As +9%) \JW+9l, 

V = \J 9s + 9% {cos ip sin 7 + sin ip cos 7 } 


/ 9$ . ! 9c 

cos ip = j = r , sin ip = 


cos 7 


v^+if 


\JW+Q 


(6 N 

ip = arc [sin 1 / 1 , cos i/>] or ip = tan -1 [ -p 


y = \l9 2 s+9l sin (7 ± ip) 

The magnitude constraint is just a constraint on the amplitude; hence the constraints become: 


6 b-i) s 

+ 

e b-i) c 

< 

A N - 1 

e N s 

+ 

d N c 

< 

a n 

6 In+i) s 

+ 

e b+i) c 

< 

a n + 1 


2.2 Nonrotating System Fourier Coefficients 

If given the Fourier Coefficients for an entity expressed in the Nonrotating System, say 9col c > 
9cOL s , 9 lat 9 lat , 9 long c i 9 lat s 


3 



where 

9 

COL 

LAT 

LONG 


c 

s 


either a collective, or a lateral cyclic, or a longitudinal cyclic coefficient 

refers to a collective coefficient 

refers to a lateral cyclic coefficient 

refers to a longitudinal cyclic coefficient 

refers to a cosine term 

refers to a sine term 


then for oscillations at m per rev (see ref. 2) 

0 = 6cOL + 0 LAT cos + &LONG sin 


where 

®COL = ®COLq + 6col s sinraflt + #cm c cosm ^ 

Olat = Olat 0 + 9 LATs sin mOt + 9 LATc cosmQt 

QlonG = @ LON Go + 6 LONGs sinmfif + 9 LO nG c c 

t is the time 


In a manner similar to that employed for the Rotating System Fourier Coefficients (see sec. 2.1), 
0COL> @LAT-> and 8 LONG can be expressed 

°COL = 0 CO Lo + \J Q'co L s + °CO L c ( cos 'POOL sin mOt + sin xp co l cos mftt} 

Hat = HaT 0 + \J e LAT s + e LAT c {cos xPlat sin mLtt + sin ipLAT cos mLlt} 

0LONG = 6long 0 + \f^ONGN^LONG c { c05x Plong sinmfif + sin^ iC>iVG cosmf2t} 

where 


and 


Then 


i’COLii’LAT, and ip LONG are phase angles 


cos ipcOL = 
c°s ip LAT = 
cos ip LONG = 


6col k 


\/ e COL s + °COL c 

Ha t s 

\/ 6 LAT s + d LAT c 

9long s 




+ 0?< 


LONGs v LONG 


v P ool = ^c[simp C0L , cos tp C0L ] 

TpLAT = arc [sin iPlat, cos iPlat] 


sin ip COL = 
sin ip LAT = 
sin ip LONG = 

or i’COL 


QcOLr 


\l e COL s + d COL c 
°LAT C 

\J e LAT s + e LAT c 
e LONGs 




+ 0?i 


LONG a ^ °LONG 


\ e COL s 

W 


‘) 


or ''Plat — tan 1 ( \ 


Hong = Mc[sinip LONG ,cosip LO NG\ or Hong = tan -1 
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The magnitude constraints on collective, lateral cyclic, and longitudinal cyclic are just constraints on 
their respective amplitudes; specifically: 


e COL s 

+ 

d COL c 

< 

a col 

e LAT s 

+ 

9 LAT C 

< 

a lat 

e LONG s 

+ 

e LONG c 

< 

a long 


2.3 Individual Actuator Stroke 

The magnitude constraints are just constraints on the vertical excursion of the individual actuators; 
specifically: 

z\ < A\ 

Z3 < M 


If the motion of the actuators is sinusoidal, 21,22 , and 23 are expressed: 



Z\ = 

2 Oi 

+ 

Ci sin (iVftt + <f> 1) 


Z2 = 

z o 2 

+ 

C2 sin (iVftt + fa) 


z 3 = 

z O3 

+ 

C3 sin(iVftf + 03 ) 

where 





N 

is the number of blades 




ft 

is the rotor rotation rate 




t 

is the time 




4>i 

is phase angle for the i-th actuator; i 

= 1 , 

2,3 

ZQ 

is the reference actuator length 





The magnitude constraints for the individual actuators becomes 


Ci < A\ 
C2 < A2 
O3 < A3 


2.4 Blade Pitch Expressed as Function of Blade Azimuth and Actuator Stroke 

To obtain an expression for the blade pitch as a function of blade azimuth and actuator stroke, the 
motion of the swashplate expressed as a function of the individual actuator strokes is derived first (see 
subsection 2 . 4 . 1 ). Next, the vertical position of a general field point (i.e., the end of a pitch link) is 
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defined (see subsection 2.4.2). Then finally, the blade pitch angle is defined from the vertical position 
of the pitch link (see subsection 2.4.23). 

2.4.1 Canonical Swashplate Plane Equation 

A necessary condition for a general field point P(x, y, z ) to lie on the swashplate plane (see fig. 1) 
is: 

(R - Rl) * [(-^2 - -Rl) X (#3 - R\)] = 0 (1) 


R = 
Rl = 

r 2 = 

Rl = 

Let V denote the matrix defined by 


OP 

= xi 

+ 

yj 

+ 

zk 

OPi 

= x\ i 

+ 

y\3 

+ 

z\k 

OP 2 

= X 2 i 

+ 

y 2J 

+ 

z 2 k 

OP 3 

= 

+ 

V3J 

+ 

z 3 k 


V = 


x — X\ 
X2 - XI 
XZ-X\ 


y-y i 
2/2 -y\ 

2/3 — 2/1 


Z — Z\ 

22 - 21 
23-21 


The swashplate plane defined by equation (1) can be expressed 


or still 


Det(P) = 


x — X\ 

— x \ 

X^-Xi 


2 / - 2/1 
2/2 - 2/1 
2/3 — 2/1 


z — Z\ 

22 -21 

23 - 21 


= 0 


C x {x - xi) + C y {y - yi) + C z {z - z\) = 0 

where 

C x = cof v {x-x 1) = 0/2 “2/l)(23 - zi) - (2/3 - 2/l)(22 “ 21) 

C y = cofx>(y-yi) = (23 - xi)(z 2 - z{) - (x 2 - xi)(z 3 - z{) 

Cz = cof v (z-zi ) = (*2 -*1X2/3 -1/1) - (X3 — Xi)(y2 — 2/1) 

Equation (4) can be rewritten as: 


Let 

and 


Then 


C x x + Cyy + CzZ — C x x 1 + Cyy\ + CzZ\ 


D — C x x 1 + Cyy\ + C z z\ 



Dy 



D z = 


D_ 

C 2 


x y z 

D x Dy D z 


Canonical form 
of the swashplate 
plane equation 


( 2 ) 


(3) 

(4) 


(5) 
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Points Pi(xi, y-j, z-j), P2(X2, Y2 , ^ 2 ), and P3(x 3 , y 3 , z 3 ) 
define the swashplate plane. 

Coordinates z-j, z 2 , z 3 define the vertical position of 
the actuator reference points on the swashplate plane 

The x-y plane is the Nonrotating System datum plane. 

The z-axis is the Shaft Axis. 

Figure 1. Swashplate plane geometry. 
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To evaluate D, note that 


D = 


Let 71 denote the matrix defined by 


Xl 

2/1 

21 


Xl 

2/1 

21 

x 2 - Xl 

2/2 “2/1 

N 

to 

1 

N 
►— • 

= 

x 2 

2/2 

22 

XZ - Xl 

2/3 — 2/1 

23 ~ 21 


x 3 

y 3 

23 


n = 


x \ 

x 2 

x 3 


y i 

V2 

y 3 


z 2 

z 3 J 


Then 


x i y i 

x 2 2/2 
^3 2/3 


z\ 

z 2 

z 3 


( 6 ) 


D = Det(1Z) = 

2.4.2 Vertical Position of a General Field Point 

The vertical position of a general field point on the swashplate plane is defined by solving equa- 


tion (5) for 2 ; specifically: 

z = Dz ( 1 -^- JL ) 

V Dx Dy ) 

(7) 

Recalling that 

Dx = 77“ J D y — — , D z = — 


Equation (7) can be expressed 

z = ^-( p - C x x - C y y) 

(8) 

Recalling that 



Cx 

= (2/2 - 2/i) (23 - 21) - ( 2/3 - 2/1) (22 - 21) 


Cy 

= (^3 “ ^l)(22 “ 21) - (x 2 - xi)(z 3 - zi ) 



Cz = ( x 2 ~ x l ) (2/3 “ 2/1 ) - (Z3 “ Zl)(y2 - 2/l) 


and 


Then 


D = Det(Jt) = 


a?i yi 

^2 2/2 22 

S3 2/3 23 
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c x = (2/3 - 2/2)21 + (2/1 - 2/3)22 + (2/2 - 2/1)23 

Cy = (x 2 - X 3 )z 1 + (x 3 - x\ )z 2 + (x x - x 2 )z 3 

c z = ( x 2 - ^1X2/3 -yi)-( x s- x i)(y2 - 2/1) 


(9) 

GO) 

(ii) 



and 


D = (z 2 2/3 - ®32/2) 2 l + (#32/1 - xi yz)z 2 + (^12/2 — X 2V\) Z 3 (12) 

Note that z\, z 2 , and 23 do not appear in equation (11). 

Equation ( 8 ) can now be expressed 

2 = [ X 2V3 ~ *32/2 ~ (V3 ~ Ifc)* - (*2 “ *3)y]*l 
+ [ x 3Vl ~ 2/12/3 “ (2/1 - 2/3)® “ (®3 - ^l)2/] -^2 

+[zi 2/2 ~ S 22/1 ~ ( 2/2 ~ 2 /l)z ~ (a?l ~ a: 2)2/]^3 (13) 

( x 2 ~ X 1 ) (2/3 - 2/l) - (*3 ~ *l)(2/2 “ 2/l) 

which has the form 

2 = (A\x + Biy + Ci)zi + (A 2 x + B 2 y + 62)22 + (^30; + B%y + 63)23 (14) 

where 



2.4.2.1 Vertical position of a general field point expressed in the nonrotating system. Equa- 
tion (14), which is expressed in the Nonrotating System, defines the vertical position of a general field 
point on the swashplate plane. If the actuators move sinusoidally with N/ rev frequency, the vertical 
position of the actuator reference points (i.e., zi,z 2 , and 23) can be expressed: 

z\ — ZQ-y 4* E\ sin (iVQt -1- 4 >\) 

22 = 2q 2 + E 2 sin(iVT)t + 0 2 ) ( 16 ) 

z 3 = z 0 3 + -E3 sin (NQ,t + 03) 

Substitution of equations (16) into equation (14) yields: 

2 = (A\x + Byy + 61) [20J + E\ sin(iVf 2 £ -I- 0 i)] 

+(A2X + B 2 y + 6 2 )[2 o 2 + E 2 sin (NQt + 0 2 )] (17) 

-f (A3X + B^y + 63)^03 + E$ s\n(NClt + 03)] 
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An appropriate Control 6- Vector is suggested; specifically: 



( 18 ) 


Equation ( 17 ) with its three sinusoidal terms can be expressed in a form which has only one sinusoidal 
term by using the Rule for the Linear Combination of Monofrequency Sinusoidal Terms which is 
described in appendix A. 


2 A. 2.2 Vertical position of a general field point expressed in the rotating system. Equation ( 14 ), 
which is expressed in the Nonrotating System, can be applied to the Rotating System by defining the 
x and y coordinates of reference points on the pitch links as they rotate about the rotor shaft axis. For 
convenience in demonstrating the harmonic nature of the kinematics, it is assumed that the pitch links 
describe a right circular cylinder as they rotate about the rotor shaft axis (see fig. 2 ). It is emphasized 
that, in general, this assumption need not be made and that the pitch link motion need not describe a 
right circular cylinder in order to obtain the harmonic kinematic form derived subsequently. For those 
cases in which the up/down motion of the pitch links can be adequately described using small angle 
approximations, the x and y coordinates of the pitch links are expressed by: 

x = rcos(ft£ + <f) s ) 

y = r sin(ft£ + 0 S ) (19) 

where 

r is the radius of the right circular cylinder 

t is the time 

ft is the rotor rotation rate 

0s is the phase angle 

If the actuators move sinusoidally with Nl rev frequency, the vertical position of the actuator reference 
points (i.e., zi,Z2, and 23) can be expressed: 

z l = 2 0i + Ei sin (JVftf + 0 i ) 

z 2 = z 0 2 + e 2 sin (NQt + 02 ) ( 20 ) 

23 = z Q z + £3 sin (NQt + 03) 


Substitution of equations ( 19 ) and (20) into equation ( 14 ) yields: 

2 = [Aircos(ftf 4 - 0 S ) 4 - B\r sin(fti + 0 S ) + Ci][z 0l + E\ sin (NClt + 0 X )] 

+[A 2 r cos(ftf + 0 S ) + B 2 rsm(£U + 0 S ) + C 2 )[zq 2 + E 2 sin (iVftf + 0 2 )] ( 21 ) 

+[A 3 rcos(ft£ + 0 S ) + B 3 rsin(fti + 0 5 ) + C 3 ][2 03 + £3 sin (JVftt + 0 3 )] 
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Figure 2. Actuator, swashplate plane, and pitch-link geometry. 
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As in the case of the Nonrotating System (see sec. 2.4.2. 1), an appropriate Control 6- Vector is: 

Ei 

E o 

E 3 
01 
02 
03 

Expanding equation (21) by taking the indicated products yields: 

2 = Ai 2 0l r cos(ftf + fa ) + Biz 0 l rsm(nt + <f> a ) + C\z Ql 

+AiE\r sin(M2f + 0 i) cos(f2f + 0 S ) + B\E\r sin(M2£ + fa) sin(ftf + fa) 

~\~C\ E\ sin(NQ,t -f- 0j) 

+^ 2 2 02 r cos(£lt + ) + B 2 ZQ 2 rsm(Q.t + fa) + C 2 Zq 2 

+A 2 E 2 r sin (Ntlt + fa) cos (tit + fa) + B 2 E 2 r sin (NQt + fa) sin(f2f + fa) > (23) 

+C 2 E 2 sin(iVf)£ + fa) 

+A 3 z 03 rcos(nt + fa) + B 3 2o 3 rsin(f2f + fa) + C 3 zo 3 

+A 3 £ 3 r sin(NQt + 0 3 ) cos(fif + fa) + £ 3 £ 3 r sin(iVQf + fa) sin(ft£ + fa) 

+C 3 E 3 sin(iVfi£ + fa) 

Noting that: 

sin(a±/3) = sin a cos /? ± cos a sin f3 

cos (a±(3) = cos a cos (3 zp sin a sin /? 
sin a sin (3 = cos(a + /?) + \ cos(a - 0) 

sin a cos /? = ^ sin(a + /?) + ^ sin(a - /?) 

cos a sin 0 = | sin(a + /3) - \ sin(a - /?) 

cos a cos/? = \ cos(a + 0) + \ cos(a - 0) 
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Equation (23) is further expanded 


2 = ^i-20i r [ cos 0s cos Qf ~ si n 0s sin Qt] + £i20i r [ cos 0s sin Qt + sin <f> s cos Qt] 
+C\zq 1 + AiJ^irfcos 0i sin iVQt + sin 0i cos iVQt] [cos 0 S cos Qt — sin 0 S sin fit] 
+£l£ir[cos 0 i sin iVQt + sin <j>\ cos iVQt] [cos 0 S sin Qt + sin 0 5 cos fit] 

+Cii?i [cos 0i sin NQt 4- sin 0 i cos NQt] 

+A2202 7, [ cos 0s cos tot ~ sin 0s sin fif] + B 2 ZQ 2 r[cos 0s sin ^ + sin 0s cos fit] 
+C2202 + ^2 j ^ 2 7 '[cos 02 sin JVfit + sin 02 cos NQt] [cos 0 S cos Qt — sin 0 S sin fit] 
+B 2 E 2 r[c os 02 sin iVQt + sin 0 2 cos TV fit] [cos 0 S sin Qt + sin 0 5 cos Qt] 

+C2 jE ?2 [cos 02 sin iVQt + sin 0 2 cos iVQt] 

+^32/o 3 r[cos 0s cos Qt — sin 0 5 sin Qt] + £32037- [cos 0 S sin Qt + sin 0 S cos Qt] 
+C3203 + A^E^r [ cos 03 sin NQt + sin 0 3 cos iVQt] [cos 0 S cos Qt — sin 0 S sin Qt] 
+£3£3r[cos 03 sin iVQt + sin 0 3 cos iVQt] [cos 0 S sin Qt + sin 0 S cos Qt] 
+C 3 £ 3 [cos 03 sin iVQt + sin 0 3 cos iVQt] 

Expanding the products for fc =1, 2, 3 

[cos 0£ sin iVQt + sin 0/- cos iVQt] [cos 0 S cos Qt — sin 0 S sin Qt] 

= cos 0£ cos 0 5 sin iVQt cos Qt + sin 0*. cos 0 S cos iVQt cos Qt 
— cos 0 sin 0s sin iVQt sin Qt — sin 0*. sin 0 S cos NQt sin Qt 
= ^ cos 0£ cos 0s {sin[(iV + l)Qt] + sin[(iV — l)Qt]} 

+2 sin 0jt cos 0s {cos[(iV + l)Qt] + cos[(iV — l)Qt]} 

+2 cos 0£ sin 0s {cos[(iV + l)Qt] — cos[(iV — l)Qt]} 

— 2 sin0jt sin 0s {sin[(iV + l)Qt] - sin[(iV — l)Qt]} 

= 5 [cos <pk cos 0s — sin 0*. sin 0 S ] sin[(iV + l)Qt] 

+2 [sin 0*. cos 0s + cos 0*. sin 0 S ] cos[(iV + l)Qt] 

+ 2 [cos 0*. cos 0s + sin 0*. sin 0 S ] sin [(N — 1 ) Qt] 

+2 [sin 0*. cos 0s — cos 0*. sin 0 S ] cos[(iV — l)Qt] 

= 2 cos (0fc + 0s) sin[(iV + l)Qt] 

+2 sin(0fc + 0s) cos[(iV + l)Qt] 

+2 cos (0*. — 0s) sin[(iV - l)Qt] 

+2 sin(0j. — 0s) cos[(iV — l)Qt] 

for k = 1, 2, 3 



and likewise: 


[cos fa sin Nflt + sin fa cos Nflt] [cos fa sin fit + sin fa cos fit] 

= cos <p cos <p s sin Nflt sin fit + sin (p^ cos fa cos Nflt sin fit 
+ cos fa sin <p s sin Nflt cos fit + sin fa sin fa cos Nflt cos fit 

= ^ cos fa cos fa {- cos[(iV + + cos[(iV - l)Qi]} 

+2 sin fa cos fa {sin[(iV + 1 )flt] - sin[(7V - l)flt]} 
cos fa sin fa {sin[(iV + 1 )flt] + sin[(N - l)fii]} 

+5 sin fa sin fa {cos[(iV + l)flt] 4- cos[(iV - l)fli]} 

= ^ [sin fa cos fa + cos fa sin fa] sin[(iV + l)flt] 

~2 [cos fa cos (p s — sin fa sin0 s ] cos[(iV + 1 )flt] 

— ^[sixup^cosfa — cos fa sin0 s ] sin[(iV - 1 )flt] 

+2 [cos <pk cos fa + sin fa sin <p s ] cos [(N — 

= 2 sin (fa + <p s ) sin[(iV + l)ft£] 

— 2 cos(<£& + fa) cos[(iV + 1 )flt] 

~2 sin(fa — fa) sin[(iV — 1 )flt] 

+2 cos {fa — fa)cos[(N — 1)17^] 

for k = 1, 2, 3 


> ( 26 ) 
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Substituting equations (25) and (26) into equation (24) yields 

2 = Ai2o 1 7’[cos<?!> s cosfit — sin 0 S sin fit] 4- Si20i r [ cos 4>s sin fit + sin 0 S cos fit] 
+C\zq x 4- ^AiEir{cos(<f>i 4- 0 S ) sin[(At 4- l)f It] 4- sin(0i 4- 0 S ) cos[(AT 4- 1 )fit] 
+ cos(0i — 0 S ) sin[(A r — l)fif] 4- sin (0i — 0 S ) cos[(N - l)fit]} . 
+^BiEir{sin((j)i + 0 S ) sin[(AT + l)fi t] — cos(0i 4- 0 S ) cos[(AT 4- l)Ot] 

— sin(0i — 0 S ) sin[(At — 1 )Ot] 4- cos(0i — 0 S ) cos[(AT - l)fit]} 

+C\ E\ [cos 0 i sin NOt 4- sin <f>\ cos NOt] 

+A 2 Zo 2 r [cos 4> s cos Ot — sin <f> s sin fit] 4- B 2 zo 2 r[cos 0 S sin fit 4- sin 0 S cos fit] 
+C2Z02 "b \A 2 E 2 T {cos(<f >2 4- 0s) sin[(AT 4- l)fit] 4- sin (<fo 4- 0s) cos[(AT 4- l)fit] 
4-cos(02 — 0s) sin[(AT — l)fif] 4- sin(02 — 0s) cos[(iV — l)fit]} 
+^B 2 E 2 r{sin((f )2 + (f>s) sin[(AT 4- l)fit] — cos((fo 4- 0s ) cos[(AT + l)fit] 

— sin(02 — 4>s) sin[(AT — l)fit] + cos(02 — 0s) cos[(AT — l)fit]} 

+C 2 E 2 [cos 02 sin NOt + sin 0 2 cos iVfit] 

4-^43 zo 3 r [cos 0 S cos fit — sin 0 S sin fit] + B^ZQ 3 r[cos 0 S sin fit 4- sin 0 S cos fit] 
+C , 3- 2 0 3 + \A 3 E 3 r {cos(03 + 0s) sin[(AT 4- l)fit] 4- sin(03 + 0 S ) cos[(AT 4- l)fit] 
+ cos(03 — 0 S ) sin[(At — l)fit] + sin(03 — <j ) s ) cos[(iV — l)fit]} 
+553£ , 3r{sin(03 4- 0 S ) sin[(AT 4- l)fit] - cos(03 4- 0 S ) cos [(AT 4- l)fit] 

— sin(03 — 0s) sin[(iV — l)fit] 4- cos(03 - 0 S ) cos [(N - l)fit]} 

4-C3S3 [cos 03 sin NOt 4- sin (f> 3 cos At fit] 



Rearranging the terms in equation (27) results in: 


2 = [Cl*0i + C 2 ^ 0 2 + C32O3] 

+r[— Aizoi sin 0s + B\zq x cos 0 s - ^2^0 2 sin 0s 

+B2ZQ2 cos <j> s - A$zq 3 sin 4> s + S3203 cos <f> s ] sin Q,t 

+r[Aizo 1 cos 4> s + B\zq x sin 0 S + ^2^02 cos 0s 

+-B2202 sin 0 s + A^zo 3 cos <j> s + £3203 sin 0 S ] cos f It 

+\r[A\Ei cos( 0 i - 0 S ) - BiE 1 ! sin( 0 i - 0 S ) + A2.E2 cos(02 - 0 S ) 

-52^2 sin( 0 2 - 0 S ) + A3E3 cos (03 - 03 ) - 5 3 E 3 sin (03 - 0 S )] sin [{N - l)fit] 
+\r[AEi sin( 0 i - <f> s ) - B\Ei cos( 0 i - 0 S ) + A2E2 sin(02 - 0 S ) 

+-B2-E2 cos( 0 2 - 0 S ) + A 3 £ 3 sin( 0 3 - <f> 8 ) - B3E3 cos ( 0 3 - 0 s )j cos[(iV - l)fli] 
+[C\Ei cos 0 i + C2E2 cos 02 + C3.E3 cos 03] sin JVfif 
+[Cii?i sin 0i + C2E2 sin 0 2 + C3JE3 sin 03] cos iVfii 

cos( 0 i + 0 s) 4 - sin( 0 i + 0 s) + A2E2 cos(02 + 0 s) 

+ 5 2 £ 2 sin(02 + 0 s) + A3.E3 cos (03 + 0 S ) + -B3E3 sin (03 + 03)] sin[(iV + l)fit] 
+\r[A\E\ sin (0i + 0 5 ) - B1-B1 cos( 0 i + 0 s) + ^£^2 sin(02 4 - 0 S ) 

-B2E2 cos( 0 2 4 - 0 s) 4 - sin( 0 3 + 0 s) - -B3-E3 cos (03 + 0 s)] sin[(iV + l)fit] , 
Note: z represents an arbitrary vertical reference position along a pitch link 


(28) 


2.4.23 Blade pitch. A simplified model of the pitch link, pitch link connecting arm, and blade 
angle geometry is presented below to illustrate the nature of the geometry. This geometry is rotor head 
peculiar and, correspondingly, each case requires individual analysis. 


It is assumed that the pitch link is parallel to the rotor shaft about which it rotates and that it 
describes a right circular cylinder as it rotates about the rotor shaft (see fig. 2). It is further assumed 
that the pitch link connecting arm is nearly perpendicular to the pitch link and to the blade pitch axis 
(see fig. 3). Finally, it is assumed that the pitch angle 6 is small enough so that 

sin# « 6 
cos 6 « 1 


Then 


1(0 + 0o) = z- zoo 


e= z_-zoo_ eo 

L 

where 

6 is the incremental blade pitch angle measured from do 
6 q is an arbitrary baseline trim value of 9 


( 29 ) 
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1 is the pitch link connecting arm length 

2 is the vertical position of the end of the pitch link 
zqq is the reference value of z which corresponds to 9q 

2.5 Rate Control 

Rate control is achieved by first taking the time derivative of the entity whose rate is to be 
constrained, and then by defining a constraint function for this time derivative. The entities whose rates 
are to be determined have one of two basic forms; these are: 

y = Asinwf + Bcosut 


and 


y = yo + c sin(ujt + 7 ) 


The first form is used for Fourier equations expressed in both the Rotating and Nonrotating Systems as 
described in sections 2.1 and 2.2, respectively. The second form is used for the sinusoidal individual 
actuator strokes described in section 2.3. The time derivative and constraint function for each of these 
forms are presented in subsections 2.5.1 and 2.5.2, respectively. 


2.5.1 Rates for the Form y = A sin ut 4- B cos ut 


If given 


Then 

Noting that 


y = A sin ut + B cos ut (30) 

y = uA cos ut — uB sin ut 


sin(a ± /?) = sin a cos /3 ± cos a sin f3 
cos(q ± (3) = cos a cos (3 ^ sin a sin /? 

y — \j M) 2 + (uB) 2 | ( — = 1 co swt - ( ) sin 


y = wVA 2 + £ 2 {cos ut cos 7 — sin ut sin 7 } 
where 7 is the phase angle defined by 

A 


cos 7 = 


sin 7 = 


B 


VA 2 + B 2 ’ ““ ' VA 2 + B 2 

7 = arc(sin 7 , cos 7 ) or 7 = tan -1 
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So: 


(31) 


y = u>y/~$ + B 2 cos (ut 4- 7 ) 

The rate constraint is just a constraint on the amplitude of equation (31); specifically: 

u\[a 2 ~+B^<A r 

2.5.2 Rates for the Form y = yo + Csin(u>f + 7 ) 

If given 


y = yo + C sin (ut + 7 ) 


Then 


y = uC cos(u )t + 7 ) 

The rate constraint is just a constraint on the amplitude of equation (5); specifically: 

uC < Aji 


(32) 


(33) 

(34) 

(35) 


2.6 Power 

The hydraulic power required to operate the actuators can be the limiting factor in the use of Higher 
Harmonic Control to reduce helicopter vibration. Correspondingly, it is desirable to be able to define a 
constraint function which limits the hydraulic power and which is expressed in terms of variables which 
are computed during the rotor simulation. Accordingly, the following procedure was defined to express 
hydraulic power in terms of mechanical variables. 

The Hydraulic Power (Pjj) requirement is roughly proportional to the Mechanical Power (Pm) 
generated; specifically: 

P H = K p P M (36) 

where K p is the constant of proportionality. 

The Mechanical Power (Pm) is proportional to the time rate-of-change of the Total Mechanical 
Energy (E) of the system; specifically: 


Pm = K E J t (E) (37) 

where Ke is the constant of proportionality. 

For the systems under consideration, the Total Mechanical Energy (E) can be assumed to be 
comprised of the Kinetic Energy (T) and the Potential Energy (V) due to the presence of an effective 
spring constant, and is the sum of the individual mechanical energies of several subsystems. These 
subsystems include nonrotating system components such as the actuators themselves and rotating system 
components such as the individual blades with hinges and linkages. The Total Mechanical Energy for 
each of these subsystems is: 

E = T + V (38) 
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The Kinetic Energy (T) is proportional to ( dz/dt ) 2 in ml 2 /t 2 . 

The Potential Energy (V) for an effective spring constant K is: 



( 39 ) 


where 

F s +K(£ — zo) 

K the effective spring constant with units of j = ^ 
z 0 vertical displacement at which the effective spring force is zero 
Zqo reference vertical displacement at which V\ is defined 


Then 


V = K f (i £-zo)df = K 
Jzoo 

(■ 




V = K \ ^z z - ZQZ + (*o*00 - 


-\Z 

zoo 

in 


mV 

T 


which can be expressed as 


where: 


V = Fz 2 + Gz + H i 


in 


ml 2 

~W 


Then 


F=\k, a = -z 0 K , H = K (z 0 zoo - izoo) 

E = D {%f + F > 2 + Gz + Hin ^ 

where D is the constant of proportionality for T 

Recalling equation (31), the Mechanical Power (Pm) can be expressed as: 

2 


p M = K e - 


D I) +^+G,+H 


, ml 2 
in — o- 
t 3 


Pm = k E 
Pm = k e 


2 D 


2 D 


( dz\ ( d 2 z s 
\ dt ) \ dt 2 , 


+ 2 f 2 (l) +G 


'dzY 

,dt)_ 



+ 2 Fz + G 


'dz\ . ml 2 
dt) m ~t 3" 


(40) 

(41) 


(42) 


(43) 


(45) 

(46) 
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2.6.1 Nonrotating Subsystems 

It is assumed that the three actuators with their associated moving linkages are the three principal 
nonrotating subsystems which must be considered. The Mechanical Power Pm for each of the actuator 
subsystems is expressed by equation (46) using the appropriate 2 position defined by equation (16) in 
section 2.4.2. 1; specifically if 

z l = z 0 i + El sin(NQt 4- 0i ) 

z 2 = z 02 + E 2 sin(NQt + fa) Equation (16) in section 2.4.2. 1 

z 3 = z 0 3 + E 3 sin(NQt + 4>z) 

then for actuator k, k = 1, 2, 3 the Mechanical Power Pm is: 

P Mk = NQE k K E (2D{-N 2 Q 2 E k sin(NQt + <f> k )} 

+2F{zQ k ■+■ Efc sin(NQt 4- .)} -1- G] cos(NQt 4* djt) (47) 

for k = 1, 2, 3 

Let sc(NQt) denote either sin(NQt) or cos(NQt) and noting that 

sin(NQt 4- 4> k) = cos <j> k sin(NQt) 4- sin <f> k cos(NQt) 
cos(NQt + = c°s 4>k cos (NQt) — sin <^. sin(NQt) 

sin 2 (NQt) = ^[1 — cos(2NQt )] 

z 

sin(NQt) cos (NQt) = ^ sin(2 NQt) 

cos 2 (NQt) = ^[1 + cos(2NQt)] 
z 

Equation (47) for the fcth actuator subsystem has the form 

= B\ sin (NQt) 4- B 2 cos(NQt) + B 3 sin 2 (NQt) + P4 sin(NQt) cos (NQt) 

+B 5 cos 2 (NQt) (48) 

Pm f. = C 0 + C\ sin(NQt 4- £j) + C 2 sixi(2NQt + £2) (49) 

where B\,B 2 , P3, P4, P5, Co, Ci, C2, Cl> an d C2 are constants. The total Mechanical Power Pm for 
all three actuator subsystems is defined by summing the Mechanical Power of each of these 
subsystems. The result has the same form as that of equation (49), but where the constants are different; 
specifically 

Pm ~ E)q 4~ D\ sin(NQt -t- ifti) 4* D 2 sin(2 NQt 4* ^2) (50) 

where Do,Di,D 2 ,ipi, and 1 P 2 are constants. 

If the power required by just the three nonrotating actuator subsystems is to be considered, then 
the power constraint can be expressed as: 

|X>oI + Pi| + |Z> 2 |<[Pm] max (51) 
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2.6.2 Rotating Subsystems 


It is assumed that the N individual blades with their hinges and linkages are the N principal 
rotating subsystems which must be considered. The Mechanical Power Pj ^ for each of the individual 
blades is expressed by equation (46) using the 2 position defined by equation (28) in section 2.4.2.2 
where the appropriate phase angle 4> $i for the 1-th blade is: 

( 360° \ 

— J for i = 1,2,..., iV (52) 


where 

N number of blades in the rotor 
i blade number; i = 1 , 2, . . . , N 

In a manner similar to that employed for the nonrotating subsystems, the aforementioned procedure 
(i.e., substitution of eq. (28) of sec. 2.4.2.2 into eq. (46)) will yield the power requirement for an 
individual blade subsystem. Combining these results for each individual blade, an expression can be 
obtained for the power requirements for all the rotating subsystems. This expression is the product of 
two terms comprised of various harmonics of fit; specifically: 

P M = (TERM!)(TERM 2 ) (53) 


where: 


TERMi = Aq + Aisc(fit) + >l 2 Sc[(iV — l)fit] + A^sc(NClt) + -f l)fit] 

TERM 2 = Bisc(fit) + B 2 Sc[(N — l)fit] + B^sc(Nflt) + B^scKN + l)fit] 


and 


Pm than has the form 


sc(-) denotes either sin(-) orcos(-) 


Pm = Cisc(fit) + C 2 sc[(N — l)fit] 4 - C3Sc(Nflt ) + C^sc[{N + l)fit] 
+C , 5sc(fit)sc(fit) - 1 - C 6 Sc[(iV — l)fit]sc(fi£) 

+CVsc(iVfi£)sc(fi£) + Cgsc[(iV + l)fit]sc(fit) 

+C9Sc[(iV — l)fit]sc[iV — l)fit] + Ciosc(iVfit)sc[(iV — l)fit] 
+Cnsc[(A r + l)fif]sc[(iV — l)fit] 

+Ci2Sc(NClt)sc(NQ,t) + Ci3Sc[(AT + l)fit]sc( 7 Vfit) 

+C'l4Sc[(A 7 ' l)fit]sc[(-/V -|- l)fit] 

The following rules and identities are used to simplify equation (55) 


(54) 


(55) 
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1. sin(a±/?) = sin a cos /? ± cos c* sin /3 
cos(a ± /3) = cos a cos j3 sin a sin / 3 


(56) 



( 57 ) 


2. sin a sin/? = -2 C0S ( a + P) + 2 Cos(a - P) 
sin a cos (3 = \ sin(a + /?) + ^ sin ( a — P) 
cos a sin /? = ^ sin(a + P) — \ sin(o: — (5) 
cos a cos /? = ^ cos(a + y0) + ^ sin(a — /?) 

3. Rule for the Linear Combination of Monofrequency Sinusoidal Terms (see app. A). 


Using the above defined rules and identities, equation (55) can be rearranged to a form with no sinusoidal 
products. 

Pm = Do + Di sin[flt + cri] + D 2 sin[2(flt) + £ 72 ] 

+£>3 sin[(iV — 2 )Qt + 0-3] + D 4 sin[(A r — l)Q,t + £ 74 ] 

+£>5 sin[iV(f}f) + £75] + Dq sin[(iV + l)fl£ + erg] 

+£>7 sin[(iV + 2)Clt + £ 77 ] + £)g sin[(2iV — 2)f It + £7g] 

+£>9 sin[(2AT — l)f It + £ 79 ] + £>iq sin[2iV(f}i) + erio] 

-KDn sin[(2/V + l)f It + cru] 4* D\2 sin[(2iV + 2)f It + £ 7 ^ 2 ] 


(58) 


where 


Table 1. 


it in equation (55) 

Sinusoid in equation (55) 

Yields sinusoidal terms of the form 

Ci 

sc(Qt) 

=^sc(f It) 

C 2 

sc[(N-l)Qt] 

=>sc[(N-l)£lt] 

C 3 

sc[N(nt)] 

=$>sc[N(SU)] 

C 4 

sc[(N+l)Qt] 

=$>sc[(N+l)nt] 

c 5 

sc(Qt)sc(£lt) 

=>sc[2(fif)] and a constant 

c 6 

sc[(N- l)Q,t]sc(ftt) 

=>sc[N(nt)] and sc[(N- 2 )Qt] 

C 7 

sc[N(Qt)]sc( tit) 

=>sc[(N+l)$H] and sc[(N-\)flt] 

C 8 

sc[(N + 1 )Oi] sc(Q,t ) 

=>sc[(N+2)Q,t] and sc[N(£lt)] 

C 9 

sc[(N-l)Qt]sc[(N-l)Qt] 

=>sc[(2N-2)Qt] and a constant 

C 10 

sc[N (£lt)]sc[(N - 1 )£lt] 

=>sc[(2iV-l)fit] and sc(CU) 

C \1 

sc[(N+l)nt]sc[(N-l)nt] 

=^sc[2/V(flf)] and sc[2(fi£)] 

C\2 

sc[N(Slt)]sc[NW] 

=*»sc[2./V(ftf)] and a constant 

Ciz 

sc[(N+mt)sc[N(nt)] 

=r>sc[(2iV+l)f2f] and sc(Qt ) 

Cl4 

sc[(N+l)£lt]sc[(N+l)£lt)] 

=^sc[(2iV+2)f2t] and constant 
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The power constraint can be written as: 


P H = K p P m < [Ph]mAX 


(59) 


K P Y \ D j I ^ \ p h\max 

j=o 

3.0 RESULTS 

Six constraint functions associated with swashplate blade pitch control were defined by expressing 
them in explicit equation form. These functions constrain: 1) blade pitch Fourier Coefficients expressed 
in the Rotating System, 2) blade pitch Fourier Coefficients expressed in the Nonrotating System, 3) stroke 
of the individual actuators expressed in the Nonrotating System, 4) blade pitch expressed as a function 
of blade azimuth and actuator stroke, 5) time rate-of-change of the aforementioned parameters, and 
6) required power. In addition to the aforementioned constraints, the kinematics of swashplate blade 
pitch control referenced to the individual strokes of the actuators in the nonrotating system were derived 
using vector analysis. 
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APPENDIX A 

RULE FOR THE LINEAR COMBINATION OF MONOFREQUENCY 

SINUSOIDAL TERMS 


Consider first the sum of three monofrequency sinusoidal terms; specifically if given 
w = F\ sin(0 + fa ) + F 2 sin (0 + <fo) + F3 sin(0 + <fo) 


Noting that 
Then 


or 


Let 


Then 


or 


or still 


sin(0 + 0) = sin 8 cos 0 + cos 8 sin 0 

w = F\ sin 8 cos <p\ + F\ cos 8 sin 0i 
+F 2 sin 8 cos (f >2 + F 2 cos 6 sin 02 
+F3 sin 6 cos 03 4- F3 cos 6 sin 0 3 

w = [Fi cos 0 1 + F 2 cos $2 + F3 cos ^3] sin 8 
+[Fi sin 0 1 + F 2 sin 0 2 + F3 sin 03] cos 8 

H s = Ficos^i +F 2 COS 02 + F3COS03 
H c = F\ sin 4>\ + F2 sin 02 + F3 sin 0 3 

w = H s sin 8 + H c cos 8 


w = 4hT+h2 


H, 




sin# + 


Hr 


7^1+ Ml, 


cos 8 


w = FF+hII cos f3 sin 8 + sin 0 cos # j 
where 0 is the phase angle defined by: 


fH \ 

0 = arc{sin 0, cos 0} or 0 = tan -1 


and 


w = Jh$ + hI sin (8 + 0) 


Now consider the general case where: 


w = F\ sin(# + 0si) 4- F 2 sin(# + 0s 2 ) + • • • + Fl sin(# + 0sx,) 

+Gi cos(# 4- 0ci) 4- G 2 cos {8 4- 0C2) + • • • + Gm cos (8 4- 0ca/ ) 


(A-l) 

(A-2) 

(A-3) 

(A-4) 

(A-5) 

(A-6) 

(A-7) 

(A-8) 

(A-9) 

(A- 10) 


25 



Noting that 


and 

Then 


sin(0 + 0) = sin 8 cos 0 + cos 8 sin 0 
cos (8 + 0) — cos 8 cos 0 + sin 6 sin 0 


w = F\ cos 0 Sl sin 8 + F\ sin 0 Sl cos 8 + F 2 cos 0 S2 sin 8 + F 2 sin 0 S2 cos 6 
+ + Fi cos 0 Sl sin 8 + Fi sin 4> Sl cos 8 

+(?i cos 0 Cl cos 8 — G\ sin 0 Cl sin 8 + G 2 cos 0 C 2 cos 8 — G 2 sin 0 C2 sin 0 
+ . . . + G^- cos 0 CM cos 0 sin 0 CM sin 0 


Equation (70) can be rewritten as 

w = H s sin 8 + H c cos 8 


where 

# s = Fi cos 0 Sl + F 2 cos <f>s 2 + • . . + Fl cos (f) SL 
—G\ sin (pc\ G 2 sin 0^2 • • • M ^tn 

H c = Fi sin 0 Sl + F 2 sin <p S2 + ... + F L sin 0 SL 

+Gl COS 0 C1 + G 2 COS 0 C2 + . . . + Gjtf cos 0 CM 
Equation (71) can be rewritten as 


or still 



w = FF+hI{ 

where (3 is the phase angle defined by: 

/ 3 = arc{sin/?,cos/?} 


cos (3 sin 8 + sin /3 cos 0 j 


or 





and 

w = + H 2 sin(0 + /?) 


(A-ll) 

(A- 12) 
(A- 13) 

(A- 14) 
(A- 15) 

(A- 16) 
(A- 17) 
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